A polaritonic two-component Bose-Hubbard model 
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We show that polaritons in an array of interacting micro-cavities with strong atom-photon cou- 
pling can form a two-component Bose-Hubbard model. Both polariton species are thereby protected 
against spontaneous emission as their atomic part is stored in two ground states of the atoms. The 
parameters of the effective model can be tuned via the driving strength of external lasers. We 
also describe a method to measure the number statistics in one cavity for each polariton species 
independently. 

PACS numbers: 03.67.Mn, 42.50.Dv, 73.43.Nq, 03.67.-a 



Introduction In recent years, significant progress in 
the theoretical and experimental study of quantum 
many-body phenomena has been made by employing ar- 
tificial structures that permit unprecedented experimen- 
tal control and measurement access. Early activity in 
this field took place in arrays of Josephson junctions [l[ 
and was followed by several important developments with 
ultracold atoms in optical lattices Despite their suc- 
cess, Josephson junction arrays and optical lattices face 
limitations as it is challenging to access and control indi- 
vidual lattice sites, due to their small separation. 

A possibility to overcome these hurdles has very re- 
cently been suggested in arrays of coupled micro-cavities, 
where a scheme for simulating the Bose-Hubbard Hamil- 
tonian Q and models of interacting Jaynes-Cummings 
Hamiltonians Q have been proposed. The phase dia- 
grams of these models were studied in Q, where the 
existence of a glassy phase has been predicted. These 
setups, where atoms interact with the resonant modes of 
the cavities, offer further possibilities to generate effective 
many-body systems as they can be manipulated with ex- 
ternal driving lasers. One of these possibilities, effective 
spin Hamiltonians, has been studied recently [fj|. 

Here, we show that coupled high-Q cavities can host 
an effective two-component Bose-Hubbard model, 
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are the polariton energies, [4, U c and l\ c their on-site 

interactions and Jb t b, Jc,c and Jb, c their tunneling rates. 

Bose-Hubbard models of two components 0] can dis- 
play several interesting phenomena which are partly also 
known for a Luttinger liquid of low energy excitations in 
fermionic systems. Among these are spin density separa- 



tion 0, Q , spin order in the Mott regime and phase 
separation 

We consider an array of cavities and study the dy- 
namics of polaritons, combined atom photon excitations, 
in this arrangement. Since the distance between adja- 
cent cavities is considerably larger than the optical wave- 
length of the resonant mode, individual cavities can be 
addressed. Photon hopping occurs between neighboring 
cavities while the force between two polaritons occupying 
the same site is generated by a large Kerr nonlinearity 
[1^ . This force can be repulsive and attractive. Each 
cavity is interacting with an ensemble of these atoms, 
which are driven by an external laser. By varying the 
intensity of the driving laser, the parameters of the ef- 
fective model can be tuned. An experimental realization 
would require cavities that operate in the strong coupling 



regime [H [H, Q M, 0, E | ■ 

The atoms To generate a force between polaritons 
that are located in the same cavity, we fill the cavity with 
4 level atoms, see figured] The transitions between levels 
2 and 3 are coupled to the laser field and the transitions 
between levels 2-4 and 1-3 couple via dipole moments 
to the cavity resonance mode. It has been shown by 
Imamoglu and co-workers, that this atom cavity system 
can show a very large nonlinearity [3] • 

In a rotating frame with respect to Hq = u>c{w a+^) + 
12j ! =i( LL 'c (T 22 + a -'C' (7 33 + 2wcoi 4 ), the Hamiltonian of the 
atoms in the cavity reads, Hj — X)^=i( £f7 22+^ <T 33 + (^ + 

324 0-24 flt + h.c.) , 
kj)(lj\ transfers level I of atom j to level 
k of the same atom, loq is the frequency of the cavity 
mode, S, A and e are detuning parameters (see figure [TJ, 
f2 is the Rabi frequency of the driving by the laser and 
gi3 and 324 are the parameters of the dipole coupling 
of the cavity mode to the respective atomic transitions 
which are all assumed to be real. All atoms interact in 
the same way with the cavity mode and hence the only 
relevant states are Dicke type dressed states 22 1. 
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FIG. 1: In each cavity atoms couple to the cavity resonance 
mode and an external driving laser. The figure shows the level 
structure and transitions of one atom, u>c is the frequency of 
the cavity mode, SI is the Rabi frequency of the driving by the 
laser, gi3 and <?24 are the parameters of the respective dipole 
couplings to the cavity mode and 8, A and e are detunings. 



Polaritons In the case where (724 = and e = , level 
4 of the atoms decouples from the dressed state excitation 
manifolds [l8| . If we furthermore assume that the num- 
ber of atoms is large, N 3> 1, the Hamiltonian Hj can be 
written in terms of three polariton species. Their creation 

(and annihilation) operators read, p = -g ( gS\ 2 — fla) 



and p ± 
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s 2 , s\ 2 = 



7^E£Li°21 and s is = ^EiLi°3i- The operators 



Pq, and pi describe polaritons, quasi particles formed 
by combinations of atom and photon excitations. 

In the relevant Hilbert space spanned by symmetric 
Dicke type dressed states and for W > 1, they satisfy 
bosonic commutation relations, [pj, p{\ = and [p^pf] = 
5ji for j,l = 0,+,—, where the neglected terms are of 
order "number of polaritons" /N. p\, p\ and pl_ thus 
describe independent bosonic particles. In terms of these 
polaritons, the Hamiltonian Hi for 924 = s ~ reads, 

[#/] 924 =o, £ =o = A*o P f oPQ + M+ P+P+ + fJ>- P-P- , where the 
frequencies are given by [1q = and fi± = (5 ± A)/2. 

We will now consider the case S 3> fi,g- Here, the 
polaritons and their frequencies read, 
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up to first order in <5 _1 . There is no spontaneous emission 
from the atomic level 2 and hence to leading order, the 
polaritons p\ and pL do not experience spontaneous emis- 
sion loss. We therefore define the two polariton species 



& f = 4 (as{ 2 - nrf) ; c t = I (ns\ 



, ga f J . (3) 

In the rotating frame, the polaritons b' have an en- 
ergy fib — and the polaritons ct have an energy 



fi c = —{B/S)B. A possible disorder in the resonance 
frequency of the cavities and hence in 8 would thus affect 
fib and fi c differently which can have interesting conse- 
quences for the phase transitions of the model [l9(. The 
dynamics of these two species is governed by the two 
component Bose-Hubbard Hamiltonian |T]) as we shall 
see. 

Perturbations To write the full Hamiltonian Hi. in 
the polariton basis, we express the operators EjLi °22 
and at EjLi °24 m terms of , ct and p^_. We obtain, 



E; =1 < 2 « « -^4 (s^c 2 - 6 2 ) + (.g 2 - n 2 )bc) /B 2 , 
where S\ 4 = ^7= EjLi u axi ano - we ma de use of the ro- 
tating wave approximation: Since 5 ^> Q,g, couplings to 
the polaritons p , are negligible, provided that 



\g 2 4\ , \e\ , |A| < \n + - fi b \ , \fJ,+ - He 



(4) 



For max.(\gugn/B%\g2i(g 2 -Sl 2 )/B 2 \) < |A|, the 
couplings to level 4 can be treated in a perturbative way. 
If furthermore \g24gfl/ B 2 \ <C \B 2 /S\, this results in en- 
ergy shifts of ri h (rib - 1) U b , n c (n c - 1) U c and n b n c U bc , 
where rib and n c are the numbers of b' respectively c' 
polaritons. The on-site interactions for the polaritons b' 
and ct can thus be written as [23| 

U b 6+6 (6+6 — 1) + f/cc+c^c-l) + U bc tfbJc (5) 

with C/ 6 = - 5 | 4 g 2 n 2 S- 4 A- 1 , C/ c = - ff 2 4 .g 2 fi 2 S- 4 (A + 
2B 2 /6)- 1 and i7 6c = -g^Cg 2 - ft 2 ) 2 B~ 4 (A + B 2 /6)-\ 
Note that C/ fc > if A < 0, U c > if A + 2B 2 /S < 0, 
Ubc > if A + B 2 /5 < and vice versa. There can thus 
be repulsive and attractive interactions at the same time, 
e.g. for A < and |A| < B 2 /S we have U b > 0, U c < 
and Ubc < 0. In a similar way, the two photon detuning 
e leads to and additional on-site term 



b)) , 



(6) 



where the transitions &t c _|_ c t^ are suppressed if 
| £5 n/i? 2 | « \B 2 /S\. 

Polariton tunneling If the cavities are either cou- 
pled by optical fiber tapers or directly via an overlap 
of evanescent fields, photons can tunnel between neigh- 
boring cavities. This process is described by the Hamil- 
tonian a (ala^, + h.c), where a is the tunneling rate 
of the photons. We translate this term into the polari- 
ton picture and assume that the tunneling rate is much 
smaller than 8. In this regime, the tunneling does not 
induce transitions between the polaritons b* or ct and 
p\ ■ Hence the p + decouple from the polaritons fct and c' 
whose tunneling terms read, 

Jbbtf-bfr + Jccctc^, - J&c(Z>tc£, + ct&#) + h.c. , (7) 

where Jbb = &g 2 /B 2 , J cc = ail? / B 2 and Jbc = agQ/B 2 . 
If \Jbc\ \B 2 /5\, transitions between and c' are sup- 
pressed. This is the case for any ft as long as g 2 ^> a5/2. 



FIG. 2: Left: The polariton interactions U b (dashed line), 
U c (gray line) and U bc (solid line) as a function of Q/g-13. 
Right: The tunneling rates \J bb \ (dashed line), |J CC | (gray 
line) and \Jbc\ (solid line) together with |/i c — Hb\ (dotted line) 
as a function of fi/313. The parameters of the system are 
324 = £713, N = 1000, A = -513/2O, S = 2000v^Vc;i3 and 
a = 313/IO. 



Parameter range Here we give one example how the 
parameters of the effective Hamiltonian ([1]) vary as a 
function of the intensity of the driving laser 0. We 
choose the parameters of the atom cavity system to be 
.924 = 513, N = 1000, A = -513/2O, S = 2000^ 5 i 3 and 
a = 913/10. Figure [2] shows the interactions Ub, U c and 
Ubc, the tunneling rates Jbb, Jcc and Jbc and |/z c — as 
a function of fi/513. For g as SI we have \Ubc\ <SC \Ub\, \U C \ 
and Jbb as Jcc ~ Jbc- Whenever \fx c — fi b \ < \J bc \, 
polaritons get converted into polaritons and vice versa 
via the tunneling Jb c - With the present choice of a and 
5, this happens for 0.16g < SI < 1.6g. To avoid such pro- 
cesses, one either needs to choose a smaller or 5 larger, 
where both choices would require higher Q of the cavities 
to ensure sufficient lifetime. The interactions Ub, U c and 
Ubc can furthermore be adjusted by varying the detuning 
A. This can be done by generating a Stark shift to the 
atomic level 4 with an additional laser that drives the 
transition between level 4 and a further atomic level in a 
dispersive (detuned) way. 

Numerical results To confirm the validity of the ap- 
proximations involved in the above derivation, we present 
a numerical simulation of the full dynamics of polaritons 

and in three cavities that each couple to N = 1000 
atoms and compare it to the dynamics of the correspond- 
ing effective model (JTJ) . We consider initial conditions 
with exactly one polariton in cavity 1 and in cavity 2 
and exactly one polariton c< in cavity 3. Figure [3^ shows 
the numbers Nb = (rib) and N c = (n c ) of polaritons 



and and their number fluctuations Fh = (n^) 



and F c = (n^.) — N^ for the first cavity. Figure [3Jd in 
turn shows differences between the full description and 
the effective model ©, [JVjfiui-^BH, [N e ] m -[N D ]BB, 
[-F&lfuii — [-F&]bh and [F c ] fuU - [F c ] B h- The effective model 
describes the dynamics very well. 

Spontaneous emission and cavity decay Level 2 of the 
atoms is metastable and hence its decay rate negligible 
on the relevant time scales. The decay mechanisms for 
the polaritons and thus originate in the cavity de- 
cay of the photons and the very small but non-negligible 
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FIG. 3: a: N b (dotted line), iV c (dashed line), F b (solid line) 
and F c (dash-dotted line) for a full model of 3 cavities with 
£/24 = 3i3, s = 0, N = 1000, fi = fvlVflia, 5 = 10 4 3i 3 , 
A = —46313 and a = —2.2 x 10~ 3 3i3. b: Differences be- 
tween the full and the effective description, [JVi,]f u ii — [M>]bh 
(dotted line), [iVc] fu n - [A^cJbh (dashed line), [F b } {uU - [F b ] B H 
(solid line) and [.F c ]f u ii — [-FJbh (dash-dotted line) for the same 
model. 



occupations of the excited levels 3 and 4. The occupa- 
tion of level 4 is due to the coupling X)jLi(°42 a + h- c -), 
whereas the occupation of level 3 only affects the po- 
laritons and stems from the linear correction term 
— (B/S)S\ 3 in equation @. The resulting effective de- 
cay rates, Tb for the polaritons and T c for the po- 
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- 2)^1^-74 and 
74, where O is the 



Heaveside step function, k the cavity decay rate and 73 
(74) the spontaneous emission rates from levels 3 (4). For 
successfully observing the dynamics and phases of the ef- 
fective Hamiltonian |T]), the interactions Ub, U c and Ubc 
need to be much larger than Tb and T c . 

The experimentally least demanding case is the one- 
component model for the polaritons b\ for which 5 ~ g. 
Assuming 1724 = 513 the maximal achievable ratio of 
Ub/Tb is here | gi3/ \/k @(n b — 2)74. In particular the 
Mott state for the polaritons w , where nb ^ 2, can even 
be realized in bad cavities without the strong coupling 
regime. However, to observe the transition to the super- 
fluid phase, the strong coupling regime with 1713 3> 
is required for the single component model, too. 

To obtain an estimate for a model with both compo- 
nents, and c\ we consider three cases, g as S7, SI as lOg 
and SI as g/10. Note that g <C S and hence spontaneous 
emission via level 3 is strongly suppressed. Denoting £ = 
.913/ \f K li, the achievable ratios of interaction versus de- 
cay rates for g as SI are Ub/Tb as U c /T c as C/(2v / 2), while 
the cross interaction vanishes, Ubc ~ 0. For SI = IO.9 
(SI = <?/10) the achievable ratios are Ub/Tb as C/100 
(U b /Tb « C/2), Uc/Tc « C/2 (Uc/Tc as C/100), and 
Ubc/ma,x(Tb,T c ) ~ ( ([4 c /max(r 6 , T c ) ~ £). 

Realizing these parameters requires cavities that op- 
erate in the strong coupling regime with large cooper- 
ativity factors, £ 3> 1. This regime is currently being 
achieved in several devices, in photonic band gap cav- 
ities [H (C ~ 3), Fabry-Perot cavities (C ~ 13), 
toroidal micro-cavities 14 1 (C at 7), fiber cavities [15[ 
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(( P 5 17) and micro-cavities on a gold coated silicon chip 
jlq (C ~ 6) among others. Our scheme should thus be 
experimentally feasible with current or soon to be avail- 
able technology. Values of £ that are predicted to be 
achievable are as high as 200 for photonic band gap cavi- 
ties and 3000 for toroidal micro-cavities 17[ . Besides the 
strong coupling, a realization of our scheme also requires 
trapping the atoms in the location of strong coupling for 
sufficient time. 

Measurements The number statistics of both polari- 
ton species and in one cavity can be measured using 
state selective resonance fluorescence in a way proposed 
in 21 1 . In the one-component BH model the polari- 



tons can therefore be mapped by a STIRAP passage [20[ 
onto the atomic levels. In the two-component case the 
STIRAP can however not be applied as in [§| because 
the energies /Z& and fi c are similar and the passage would 
thus need to be extremely slow to still be adiabatic. 

For two components, one can do the measurements as 
follows. First the external driving laser fl is switched 
off. Then the roles of atomic levels 1 and 2 are inter- 
changed in each atom via a Raman transition by apply- 
ing a 7r/2-pulse. To this end the transitions 1 <-> 3 and 
2^3 are driven with two lasers (both have the same 
Rabi frequency A) in two-photon resonance for a time 
T = tt(5a/|A| 2 ((5a is the detuning from atomic level 3). 
The configuration is shown in figure [4^,. This pulse re- 
sults in the mapping \lj) <-> \2j) for all atoms j. 

Next another laser, 9, that drives the transition 1 <-> 4 
is switched on, see figure |4)d. Together with the coupling 
<?24 j this configuration can be described in terms of three 
polaritons, %, q\_ and gl, in an analogous way to pj, p* + 
and p_ , where now the roles of the atomic levels 1 and 
2 and the levels 3 and 4 are interchanged. Hence, if wc 
choose = the ir/2- pulse maps the onto the dark 
state polaritons of the new configuration, gj, whereas if 
we choose 6 = — f2 it maps the c< onto q\. The driving 
laser is then adiabatically switched off, G — > 0, and the 
corresponding STIRAP process maps the gj completely 
onto atomic excitations of level 1. This process can now 
be fast since the detuning A is significantly smaller than 
8 and hence the energies of all polariton species gj, q\_ 
and qL well separated. Another tt/2- pulse finally maps 
the excitations of level 1 onto excitations of level 2, which 
can be measured by state selective resonance fluorescence 
in the same way as discussed in [3 , 21 1 . 
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FIG. 4: a: Configuration of the 7r/2-pulse. Two driving lasers 
in two-photon transition with identical Rabi frequencies A 
couple to the atomic transitions 1 <-> 3 and 2^3. b: Config- 
uration for the STIRAP process. A driving laser couples to 
the 1 <^ 4 transition with Rabi frequency Q. The cavity mode 
couples to transitions 2^4 and 1 <-> 3, where the coupling 
to 1 <-> 3 is ineffective and not shown. 



pled arrays of high-Q cavities. As a new feature, the 
model can display transitions between the two particle 
species. An experimental realization is feasible with cav- 
ities that have cooperativity factors much greater than 
unity and interact with the atoms for sufficient time. The 
local particle number statistics of both species can be 
measured independently with high accuracy. 
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SRC grant EP/E058256, the Alexander von Humboldt 
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Cientffico e Tecnologico (CNPq) and the Royal Society. 
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The whole sequence of 7r/2-pulse, STIRAP process and rgi 
another 7r/2-pulse can be done much faster than the 9 
timescale set by the dynamics of the Hamiltonian (JT|) Q 
and or can be mapped onto atomic excitations in a 

[10] 

time in which they are not able to move between sites. 

The procedure thus allows to measure the instantaneous , , 

local particle statistics of each species separately. 

Summary We have shown that a two-component [12] 
Bose-Hubbard model of polaritins can be created in cou- [13] 
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